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Abstract: We extend the effective field theory treatment of the thermodynamics of small 
compactified black holes to the case of charged black holes. The relevant thermodynamic 
quantities are computed to second order in the parameter A ~ {r^/ L)'^~^ . We discuss how 
the addition of charge to a caged black hole may delay the phase transition to a black 
string. In the extremal limit, we construct an exact black hole solution which serves as 
a check for our perturbative results. Finite size effects are also included through higher 
order operators in the worldline action. We calculate how the thermodynamic quantities 
are modified in the presence of these operators, and show they enter beyond order as in 
the uncharged case. Finally, we use the exact solution to constrain the Wilson coefficients 
of the finite size operators in the extremal limit. 
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1. Introduction 

The physics of black objects in higher dimensional spacetimes can be much richer than 
in four dimensions. In higher dimensions, it is possible to find black objects of different 
horizon topologies and transitions between these different phases. One extensively studied 
example can be seen in the rotating stationary solutions of [||], where both black hole (BH) 
and black ring solutions exist. For particular combinations of mass and angular momentum 
both phases can exist simultaneously and this has lead to a rich phase diagram. 

Another important system with black objects of different topologies and interesting 
transitions is the BH-black string (BS) system; see Q and references therein. Here a 
spacetime with a single compact dimension is considered, namely R'^"^'^ x S^, where d is 
the total spacetime dimension. Beyond M'^^^'^, this is the simplest possible extension to 
the spacetime manifold, and thus presents a natural candidate for study. In this paper 
we also calculate in this spacetime. Note that in the case of a more elaborate manifold 
compactification, one expects the transition physics to be similar to the case. 
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In the BH-BS system, the BH and BS each possess distinct topologies - the BH horizon 
has a spherical topology S*^"^ and the BS has a horizon topology S*^^^ x S^. The stability 
of both phases depends on the ratio of the two scales in the problem, the mass M and 
the size of the compact dimension L. For large masses GM ^ L'^^^, a uniform BS is 
the only stable solution. The analytical solution in this region of phase space is well- 
known and exhibits a Gregory-Laflamme instability at lower masses j^. For small masses 
GM <C L'^~^ however, the only stable phase is a BH. Since this BH is localized in a 
compact dimension, it is referred to ged BH or Kaluza-Klein (KK) BH. Moreover, 

in an intermediate mass region, there is an extra phase corresponding to a non-uniform 
BS, where the non-uniformity exists in the compact dimension. This phase may be stable 
or unstable depending on the total number of spacetime dimensions. 

Recently, the phase transitions between the uniform BS, the non-uniform BS, and the 
caged BH have been discussed extensively. In particular, numerical research of the BH/BS 
system has shown the BH and the non-uniform BS phases merge at a topology changing 
transition Although no closed form metric exists for a caged BH, it is possible to 
calculate the properties of small caged BHs in a perturbation expansion in the parame- 
ter A ~ {r(j/L)'^~^. Here rg is the Schwarzschild radius and L is the size of the compact 
dimension. Note that in the region where the BH/BS phase transition occurs, pertur- 
bative methods break down and numerical methods or other non-perturbative analytical 
approaches Q must be employed. 

There exist several methods to compute the perturbative expansion. In Q , this prob- 
lem was studied within adapted coordinates in a single patch. In Q, the thermodynamics 
were computed by applying the matched asymptotic expansion method. Finally ^ used 
an effective field theory (EFT) treatment to calculate the thermodynamics of neutral caged 
BHs. This last approach appears to be the most systematic and transparent and we em- 
ploy it here. It applies the worldline EFT methods used to describe extended objects in 
general relativity (GR), originally proposed by Goldberger and Rothstein [^, |lO| (see |11| 
for a pedagogical introduction). 

In the EFT method was used to study the thermodynamics of static caged BHs 
without charge up to ©(A^) in the perturbation expansion. The method was further opti- 
mized in where a metric parametrization using a temporal KK reduction was intro- 
duced. This parametrization improves EFT calculations in static or weakly time dependent 
systems [13|. The leading contributions to the thermodynamics of small rotating caged BHs 
were also studied in fl^ . 

In addition to the thermodynamic calculations mentioned above, the worldline EFT 
approach to GR has been successfully applied to the post-Newtonian expansion of binary 
systems P, 13|. Spin degrees of freedom have also been incorporated into the EFT frame- 



work 1 14, and gravitational radiation has been studied [^, [T^. The self force on a 
compact object was also computed and absorptive effects were introduced in |18, 19|. 
Finally, the energy momentum tensor was considered in |20| and recently, |21] showed how 
to constrain the three and four graviton vertices. 

In this paper, we extend the perturbative studies of and |12] to electrically U{1) 
charged caged BHs. Little is known about the phase diagram of this system, but a richer 
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phase structure can be expected, due to the presence of the charge as an additional parame- 



ter. The stabihty of the uniform charged BS has been considered |22] and similar properties 



to the uncharged case were found. Recently, |23] used a Harrison transformation to study 
charged black objects on KK spacetimes in Einstein-Maxwell-dilaton gravity. 

This work contains three main parts. First we compute all thermodynamic quantities 
to ©(A^) in the presence of charge. Our results are valid for any value of the charge between 
zero and its extremal value as long as A is small. Note that to O(A^), the corrections stem 
purely from the point particle description of the BH in the EFT framework. For vanishing 
charge we recover the results of [§, 12 1. 



Next, we construct an exact charged caged BH solution in the extremal limit and 
compute its thermodynamics. In this limit, the mass, tension, and electrostatic potential 
do not get renormalized, i.e. they are simply given by their expressions in uncompactified 
space. These results provide a further check of our perturbative calculations. 

Beyond 0{X^), finite size effects are expected to become important and we dedicate 
the last section of our paper to them. In the EFT approach, finite size effects are simply 
represented by higher order terms in the worldline action. We evaluate the contributions 
of the leading order finite size operators to various thermodynamic quantities for both 
uncharged and charged caged BHs. The final results are expressed in terms of undetermined 
matching coefficients, which need to be computed via a matching procedure. Of particular 
note is the uncharged caged BH, where there is a non-renormalization of the entropy at 
leading order in the finite size operators. Finally, for a small charged caged BH in the 
extremal limit, we undertake a matching calculation of the finite size coefficients using our 
exact solution. Due to a degeneracy in the thermodynamic relations, we are unable to 
uniquely fix all three coefficients in the extremal limit. 

2. Effective field theory setup 

In this section we will establish our conventions and the main ingredients for our calculation. 
The general prescription of how to set up an EFT for small KK BHs is explained in and 



12] and we refer the reader to these papers for the details such as the power counting. 

We restrict our attention to the simplest compactification in higher dimensional space- 
time, that of R'^^^'^ X S^. In this spacetime a hierarchy of scales forms when we consider a 
small BH with Schwarzschild radius tq <^ L, where L is the asymptotic size of the compact- 
ified dimension. With such a hierarchy in place, caged BHs can be treated within an EFT 
that allows their properties to be calculated perturbatively in the expansion parameter A, 
which is defined by 



^■■=[f) C(d-3), (2.1) 
with C,{z) denoting the Riemann zeta function. The Schwarzschild radius tq is given by 

._3 _ leTrGmp _ 2^V 

where mg is the mass of the BH in an uncompactified d-dimensional spacetime. In the 
language of quantum field theory mo is called the bare mass of the BH. The charge of the 
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BH is Q. Coordinates on the d-dimensional spacetime are = (x'^,x) where x = (x_|_,2;), 
with z labeling the coordinate along S^. We let Greek indices run over all coordinates while 
Latin indices denote spatial components. 

The setup of the EFT now proceeds as in ^] but we must also include the U{1) 
gauge interaction. Instead of obtaining a BH solution by solving the vacuum Einstein- 
Maxwell equations, we integrate out the BH, i.e. we integrate out the short distance scale 
tq. This is achieved by replacing the BH by an effective theory consisting of the worldline 
of a massive charged particle coupled to the Einstein-Maxwell system. The action is then 
given by 5 = ^eh-em + -^bh with 

5eh-em[5/..,«/.] = j (fx^Mim + rff^u) , (2.3) 

SB}i[x,gf_tu,a^] = -mo j dr + Q j dx^a^ + ... , (2.4) 

where f^y is the Maxwell field strength tensor, is the electromagnetic vector potential, 
and dr = gij,udxi^dx'^ . While the first two terms in 5bh describe a charged massive point 
particle, the ellipsis denote higher order terms which encode finite size effects. These terms 
are discussed in Section ^. 

Following [0] we exploit time translation symmetry and apply the standard KK ansatz 
for the metric 

ds'^ = e^'^idt - Aidx'f - e-'^^'^'^-^^-i.jdx'dx^, (2.5) 

reducing over the non-compact temporal coordinate. We can write ^ij = 6ij+aij so that flat 
Minkowski spacetime is given by the ground state (p = Ai = aij = 0. This parametrization 
reorganizes the dynamical components of the metric, i.e. g^i, {(f), Ai,aij), into a form 
which is advantageous from a computational perspective ||l2|, In addition, we also 
decompose the electromagnetic d-vector potential into temporal and spatial components 
= {ip,ai). Note the distinction between the notation for the KK gravitational scalar 
field (p and the electromagnetic scalar field ip. 

It is convenient to parameterize the worldline using coordinate time t = x^. Since our 
problem is static, we use the d- velocity dx^^/dt = (1,0) and we may fix the location of 
the BH at X = without loss of generality. As a consequence, the fields Ai, Oj, and aij 
cannot couple to the worldline as long as we neglect finite size effects, and we will drop the 
dependence on x in Sbh- Staticity also implies there is the additional symmetry of time 
reversal symmetry T, which means both vector fields Ai and a^, which are odd under T, 
can be set to zero in any purely classical calculation like the one here. In the static limit 
the action is proportional to J dt and we will suppress this component from this point. 
The resulting action takes the form 

SEH-EM[<^,a,,-,v'] = d^'-'^^\(^-R[^] + ^{d^f-2e-^^{d^f^ , (2.6) 

SBYi[^,^] = -m^e^ + Q^ + ... , (2.7) 

where {dcf))^ = {dicp){dj(f>)j^^ and similarly for (df)'^. We do not have to gauge fix the 
action because the only propagators used in our calculations are the (j) and ip propagators, 
and these can be easily computed from (^ 
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Figure 1: Feynman rules needed for the charged caged BH thermodynamics computations to 
O(A^). The single dashed lines correspond to the KK scalar field (/>, which is the propagating scalar 
mode in the gravitational sector and the double dashed lines denote the electromagnetic scalar i^. 
The Feynman rule (a) is given in momentum space, where we suppress the delta function, and 
(b)-(d) are in coordinate space. 



The action we introduced above describes the dynamics at length scales larger than 
tq. In order to compute an observable measured at infinity, we now have to integrate out 
the scale L. This is performed by splitting up the fields into short and long wavelength 
modes, e.g. we write (p = (p + (pi where (j) is the long wavelength mode and i^l is the short 
wavelength mode which we will integrate out. We will drop the index L from now on. The 
formal expression for the effective action is given by 



and clearly Ses will contain the standard Einstein-Maxwell action SEH-EMlfp, The 
BH dynamics is contained in 5efr,BH) which in practice is expanded in powers of the long 
wavelength modes and the terms in this expansion are computed via Feynman diagrams. 
These terms are then related to the required observable. For example, the mass is simply 
given by the constant term in 5ofr,BH without any couplings to the long wavelength modes. 
In [12 1 it was shown the KK parametrization of the metric in ( |2.5D reduces the number 



of diagrams to be calculated in the uncharged case. This remains true when we consider 
charged compactified BHs, however there will be diagrams at 0{X'^) which are not of an 
explicitly factorized form. These diagrams arise due to a vertex with two electromagnetic 
scalar fields ip and one gravitational field (p derived from (|27 



We now proceed to write down the Feynman rules derived from the action (2.6) and 
( [2.7] ). The long wavelength modes are not needed until Section ^ so we set them to zero for 
now. The procedure for deriving the Feynman rules is discussed in js], |l2| and we simply 
summarize the results in Fig. |^. Note we follow the conventions of |12] when deriving the 
Feynman rules. In contrast to the gravitational couplings, it is only possible to have one 
photon coupled to the worldline at a given interaction point and we also have a 4"^^ vertex 
as mentioned earlier. The propagator in position space for the gravitational (p field is 



8TTGd-3 ^ f d'^-^k 



Z)^(x^-x^;z-z) = — — — — 2^ 



2k_L-(x— x')_L+27rm(2— 2:')/L 



L d-2^^^J (2^)'^-2 k2_ + (27rn/L)2 

(2.9) 
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and for the electromagnetic scalar ip the propagator is 
D^iyL± - x.'j_; z - z') = — ^ / 



AttG sr^ f d'^ '^ki 1 gjkx-(x-x')x+27rm(2-2')/-f' 



n=— oo 



(27r)'^-2 + (27rn/L)2 

(2.10) 

In the Feynman diagrams, single dashed hnes correspond to the gravitational scalar field 
4>, whereas double dashed lines denote the electromagnetic scalar (p. The worldlines, which 
are denoted by solid double lines, do not propagate so there are no propagators associated 
with these lines. 

3. Thermodynamics of charged caged BHs from EFT 

To compute the thermodynamics of small charged caged BHs we first calculate diagram- 
matically its ADM mass and the redshift factor. The mass, for instance, is computed 
by integrating out short wavelength modes of order L, as described in (^.g|). Once these 
quantities are at hand, one can calculate the remaining thermodynamics. 

3.1 The ADM mass at 0{X^) 

The ADM mass m can be read off from the constant term in the effective action 5'effBH 



12] and to 0{X ) this equals the sum of classical vacuum diagrams shown in Fig. y. In 



the uncharged case only diagrams in Fig. |^(a) and {b) contribute and these diagrams are 



calculated to be ||12| 



Fig. 1(a) =+^A, (3.1) 



mo 2 

where the parameter A is given by (|2.1| ) and we have used the Feynman rules of Section |2|, 



Fig.|(6) = -^A^ (3.2) 



along with ( |A.1| ). 

For the charged BH there are two additional diagrams to compute. Fig. ||(c) and (d). 
The computation of diagram Fig. ^(c) is similar to Fig. |2|(a) and the result is 

It can be understood in terms of the total electrostatic energy E = ^^iQ y^i{0), where 
the sum runs over all images of the BH in the covering space and fi{0) represents the 
electrostatic potential created by the image i at the location of the BH. 

Diagram Fig. ^(d) is slightly more complicated, so we calculate it here explicitly. The 
symmetry factor we divide by is 2. We can now write the amplitude 



,3.4) 



d — 3 y—f 2 P]^ • P2 

where p = (p_l, and 



1 s-^ f d^ P-L 



L ^ J (27r)'^-2- 
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Figure 2: Diagrams that contribute to the ADM mass to 0{X^). 



By using the identity 2 • p2 = (pi + P2)^ ~ Pi ~ Pi ^^'^ then subsequently redefining the 
momentum, we find this amphtude can be written in a factorized form 



+(87rGQ)- 



d-3 



■ mo 



d-2 

and with the identity ([A.l|) it leads to 



1 

2L 



E 



1 



(27r)°'-2 p2 + (27rn/L) 



(3.5) 



Fig. = + 



Q"^ d-2 
4mo (i — 3 



(3.6) 



Combining these results gives the following expression for the ADM mass to second 
order in A 

1 



m = mo 



(3.7) 



The variable r] in the above expression plays the role of the order parameter and is given 

by 



L ld-2Q^ 
^ y 2d- 3ml' 



It vanishes in the extremal limit when 

1 d- 2 



2d-3 



Q 



2 2 

- mo 



(3.8) 



(3.9) 



holds and reaches its maximum value, rj = 1, for an uncharged BH. Note also, corrections 
to the mass of an extremally charged BH vanish to second order identically, namely the 
ADM mass is left unrenormalized to 0{X'^) in the extremal limit. As we show in Section]^, 
it turns out such a non-renormalization of the ADM mass must hold to any order in the 
perturbation expansion. 

3.2 Thermodynamics 

In the previous subsection we computed the BH mass ( p.?! ) in terms of the charge Q, 
the compactification parameter L, and the bare mass mo, however this relation does not 
fully determine the thermodynamics of the system. In this section we will compute the 
remaining thermodynamic quantities, the temperature T, the entropy S, the electrostatic 



potential and the tension r [24|. We also give the expression for the Gibbs potential G. 
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(a) (6) 
Figure 3: Diagrams tliat contribute to the asymptotic value of the electric charge. 



First, we will clarify the notion of the charge Q. In our calculations there is no 
distinction between the bare and asymptotic values of the electric charge Q. Intuitively such 
a non-renormalization stems from the electric flux conservation requirement. Technically, 
the asymptotic value of the electric charge can be read off the following linear term in the 
effective action 

5'efr,BH[0, ^ij, <f] 13 +Q(f . (3.10) 

This term is the sum of Feynman diagrams like those of Fig. ^ with (p set to a constant, 
in which case the momentum flowing into the diagram vanishes. All couplings involving 
gravitons and photons are derived from the second term in ( |2.3[ ) . Since all of these vertices 
contain derivatives acting on the electromagnetic scalar field (p, they all vanish in the 
zero momentum limit. The only possible vertex for a photon coupling to the worldline 
is Fig. |^(b) which cannot include any additional gravitons. Thus only the leading order 
diagram is nonzero and the sub-leading diagrams such as Fig. |3|(b) vanish. This means 
the asymptotic value of the electric charge is equal to its bare value Q. The inclusion 
of finite size effects does not alter this conclusion, since all higher order operators with 
electromagnetic interactions include derivatives acting on (p. 

Next we compute the asymptotic temperature. It is canonically conjugate to time 
and hence transforms inversely to t, namely T/Tq = to/t. Therefore one can relate the 
asymptotic temperature T to the local temperature Tq using the redshift. 



RTo. 



(3.11) 



As long as we do not include finite size effects, the horizon will not be deformed and the 
local temperature will coincide with the temperature of an uncompactified d-dimensional 
Reissner-Nordstrom BH |2l 



Tb(mo, Q) = T{L 



oo 



d — 3 rriQ 
d-2~S^ 



V: 



(3.12) 



where rj is the order parameter defined in (3^) and Sq is the entropy of an uncompactified 
d-dimensional charged BH given by 



So{mo,Q) 



^ d-2 



„d-3 



SirGrrii 







■(1 + ^)- 



(3.13) 



+ {d- 2)na^2 

Note the local temperature Tq in ( 3.12| ) and the expression for 5*0 in ( 3.13| ) are correct away 
from the extremal limit where subtleties might arise in the definition of temperature and 
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Figure 4: Diagrams required to compute tlie redshift factor R to 0{X^). Here, 0(0) denotes that 
the (j) field propagates back to x = on the BH worldline but does not have a vertex. 



entropy |26|. The redshift factor R in ( 3.11| ) is given by 

to 
t 



R 



(3.14) 



where (?!)(0) is defined in terms of the Feynman diagrams given in Fig. ^. We now calculate 
these diagrams to obtain R. 

The first two diagrams, Fig. ^(a) and (6) are charge independent and their value is 
given by [12 1 



Fig. ia) 
Fig. 1(6) 



(3.15) 
(3.16) 



Evidently, the computation of the redshift factor is only modified at ©(A^) by diagram 
Fig. §(c). The steps to compute this additional diagram are similar to those used to 
compute Fig. §(d). It is merely obtained by dividing the result of Fig. ^{d) by —mo- This 
in turn yields 

Q2 d-2. 



Fig. |(c 



Combining the above results, we obtain the redshift factor 

^2- 



i? = 1 - A+ 1 + 



A^ + ... . 



(3.17) 



(3.18) 



In the extremal limit r/ = 0, the redshift factor agrees with the exact solution constructed 
in Section ^. For the temperature to 0{}?) in the case of a charged compactified BH the 
result is then 



To 



1-A+ 1+ 



(3.19) 



One can close the thermodynamic equations by considering the entropy of the charged 
caged BH. Since the contribution of finite size operators are always beyond 0{}?) in the 
perturbative expansion, we can deduce horizon deformations due to the presence of the 
compact dimension can be neglected. In other words, the BH horizon is unaltered by the 
compactified dimension to 0{}?) and its area is fixed by the values of mo and Q only. Since 
the entropy of the BH is proportional to the area of the horizon, we can write down the 
entropy of the charged caged BH to 0{X^) as follows 

5o(1 + OA + OA2 + ...), (3.20) 



S 



-9- 



where So is the entropy of an uncompactified d-dimensional charged BH given by ( 3.13| ). 
In Appendix ^ we present an independent derivation of the thermodynamic relations using 
the Helmholtz free energy, which does not use the uncompactified entropy as an input. 
From the first law of thermodynamics 

dm = TdS + ^dQ + fdL, (3.21) 

one can derive the expression for the electrostatic potential 

^_,'dm\ _ d{m,S,L) (d{Q,S,L) 



dQ J S I d{Q,mo,L) \d{Q,mQ,L) 



dm\ ( dm \ / dS\ ( dS 



-1 



y \dm() J \dQ ) \dmQ ' ' (3.22) 

where we have used the standard change of thermodynamic variables in terms of Jacobians 
1 27]. Substituting explicit expressions for S and m from (3/7) and (|3.13D respectively, we 
obtain 

$ = $0 [l + ?/A-|(2-r?)A2 + ...] , (3.23) 

where ^>o = (m,o/Q)(l — rj). 

To compute the tension f, we use the Smarr relation 

{d-3){m-^Q) = {d-2)TS + fL, (3.24) 

which gives 

— = -{d- 3)r/2A -{d- Z)ifX^ + ... . (3.25) 
mo 2 

We have explicitly checked this result for the tension via a direct diagrammatic computation 
with an external cj^^, see Section ^ for further details on the procedure. 
The Gibbs free energy 

G = G{T,^,L) = m-TS -^Q, (3.26) 

which implies dG = —SdT — Qd^ + fdL. Although we did not need to use the Gibbs free 
energy above, all the thermodynamics can be derived from G by forming the appropriate 
differentials. The result reads 



G{mo,Q,L)=Go 



(3.27) 



where Go = rimQ/{d — 2) is the Gibbs potential of a charged BH in the uncompactified 
spacetime M*^"^'^. 

Examining our thermodynamic calculations, we recover the results for an uncharged 
caged BH of ^] in the vanishing charge limit r/ = 1. On the other hand, when the 
extremality condition t] = holds, we have found to ©(A^) the mass and the electrostatic 
potential remain unrenormalized and the redshift R is given by ( p.l8| ) with rj = 0. In 
addition, according to ( 3.25| ), the tension of an extremally charged BH vanishes to second 
order identically. Physically, this result states that electrostatic repulsion exactly cancels 
gravitational attraction to this order in A. These results are in agreement with the exact 
solution we construct for an extremally caged BH in the next section. 
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4. Charged caged BHs in the extremal hmit 



It has been known for some time that Einstein-Maxwell theory admits static multi-BH 
solutions [^]. The existence of these solutions can be attributed to the balance of elec- 
tromagnetic and gravitational forces when the extremal condition ( ^.Sj ) is satisfied. In this 



section we follow |2£] and and build the thermodynamics for our extremally charged 
caged BH. This will allow us to verify the perturbative expansions of the previous section 
in the extremal limit. 

4.1 Exact solution 

The action we employ is the Euclidean Einstein-Maxwell action^ supplemented with the 



Gibbons-Hawking boundary term |29|. This additional term is included to obtain an action 



depending only on the first derivatives of the metric. We therefore have 

Sm.,A,] = Y^I d'x^ (7^bJ:J - T^^J^,,) -^l^-^i^- ^o), (4.1) 

where gfi,/ and Afj, are the full metric and the full electromagnetic d- vector potential includ- 
ing all length scales involved in the problem; TZ is the Ricci scalar and J^^u = dfj_Au — d^A^ 
is the field strength tensor, hij represents an induced metric on dA4, and IC and /Co are 
the trace of the second fundamental form of the boundary for a particular metric and a 
flat metric respectively. 

For the caged BH setup, the gravitational field can be thought to originate from the 
BH and its associated images in the covering space. In other words, the setup is equivalent 
to the case in which an infinite number of identical equidistant BHs lie on a common axis 
(the z-axis) of a d-dimensional uncompactified spacetime, see Fig. |5|. This arrangement 
has the following solution |30| 

ds^ = gf^^dTf^dx" = H-^dT^ + /f^/Cd-s) (^^^2 ^ ^2^^2 ^ ^2 g-^2 Q , (4.2) 

where H extremizes the action ( [1.1| ) and is given below, dr2^_3 denotes the metric on the 
unit {d — 3)-sphere and r corresponds to the Euclidean time. The Maxwell field strength 
IS T = dA, where 

The function H is harmonic on M'^^^, with poles corresponding to the locations of the event 
horizons of the BHs and is given by 

^ " ^ + ^ J^^ (xi + {z- nLf Y^-^)/^ ' ^^'^^ 
where |x_|_| = psmO and z = pcosO. 



^In this section we change our notation to emphasize that we are not working in the EFT but in full GR, 
namely all scales are present at all times. Also, since our problem is static, switching between Lorenzian 
and Euclidean signatures presents no difficulties and is given hy t it. 
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z 



Figure 5: Schematic of exact solution construction. 



A comment should be made regarding the gauge choice of the 1-form A. As can be 
seen from ([4.3| ), = at any event horizon. Since dr is not defined on the horizon, this 
choice of gauge avoids the possibility of a singularity in the 1-form at any event horizon. 
The potential at infinity does not vanish in this gauge, so the electrostatic potential is 
given by 



$ = A(p = o)-A(oo) 



d 



-sign(Q) , 



(4.5) 



2(d-3)" 

where we have chosen to evaluate ^ at the horizon located at the origin. 

If L ^ cxD in ( |4.4| ), so that only the n = term remains, we recover the standard 
Reissner-Nordstrom extremal solution written in the isotropic radial coordinate p. Its 
relation to the Schwarzschild coordinate r is given by 

l/{d-3) 



1 



,d-3 



(4.6) 



4.2 Thermodynamics 



In what follows we evaluate the Euclidean action (4J.) on the solutions ( [4.2| ) and ( [4.31) . We 
find the Euclidean action vanishes and this leads to the conclusion that the Gibbs potential 
of an extremal charged caged BH equals zero, via the relation 



-(3G{mo,L) = \nZ 



(4.7) 



where /3 is the inverse Euclidean temperature and Z is the partition function. As shown 
in |2^, due to the topology change of the spacetime manifold, the Euclidean extremal 
Reissner-Nordstrom solution has two boundaries: at infinity (as in the non-extremal case) 
and at the horizon. We now evaluate the contribution these boundaries make to the 
Euclidean action. 

We first consider the boundary at infinity and we choose r± = const, where r± = |x_|_|, 
to represent this surface. This boundary is compact and has the topology x x E>'^~^. 
The series in (4^) can be summed explicitly if d is odd, that is d = 2j + 5, j = 0, 1, 2, . . . . 
In this case 



H = 1 + 



-IV d^ 



1 + r 



d-3 



di^lV 

-IV d^ 



E 



1 



xj_ + (z - nL)2 



2L 



d{^ly 



1 



sinh 



27r|x^| 



1 + 



2(d-4) J7d_3 L 



|x±| 



l-'^-'-l cosh 

d-4 



27r X I 



COS 



(4.8) 
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where in the last equahty we suppressed exponentially small terms of order ~ exp(— L/|x_|_|). 
These terms do not contribute further in our computations. Further, since the leading or- 
der term contains no information about parity, one can show this result is valid for any d. 
Therefore we have 

-2 

,2 



,2 / 1 , d-3 nd_2 ro f ^ 



_4X 2/{d-3) 



r=iA=-J '"-y-^' l^friV" A *. (4.9) 



One can evaluate the corresponding Gibbons-Hawking boundary term relying on the 
identity 

[ VhIC = ^ / Vh, (4.10) 
JdM c'ra Jqj^ 

where {d/dn) Jqj^ Vh is the derivative of the area of dAi as each point of the boundary is 
moved an equal distance along the outward unit normal n. Thus in our case we get 

/ Vh)c = i7-v(<i-3)^ (/5L^^,_3i/V(<i-3)^d-3^ 

J r I =const d,r±_ V / 



lr_i_=const 



(4.11) 



For fiat space /Co = (d - 3)i?-V(<i-3)^j\ therefore 

VhfCo = I3L fld-sid - 3)r'^-^. (4.12) 



/ 

Jr I 



/ r± =const 

Combining these relations together yields 

/3mo 



SttG J^^^f^^^gf d — 2 



(4.13) 



We now turn to the boundary at the horizon, given by p = 0. To compute this boundary 

contribution, we consider the hypersurface p = const and then take the limit p ^ at 
the end of the calculation. An outward-pointing normal vector for such a hypersurface is 
5^p (no sum on p) and the induced metric h^^ is given by 

and the covariant derivative of ri/^ is 

rii.; tJ, = - V^pp (^MP ^up dp In + 5^0 5^,^ dg In ^/g^ - F^^) . (4.15) 
We can now evaluate the trace of the second fundamental form on the boundary 



1 dpH 
d-3 H 



(4.16) 
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From the definition (^^J) of H, we get tlie following expansion in the vicinity of the horizon 



2 Vp 

M = ^ (l + oip^-^)) . (4.17) 
The regulator in this case is given by /Co = -{d - 2)F-V{'^-3) /p and we therefore obtain 
^ / Vh{IC-ICo) = ^ p'-' (l + . (4.18) 

OVrG Jp=const ottG V / 

As a result, the horizon boundary term vanishes when p tends to zero, that is 

^y^_^^(^-^o) = 0. (4.19) 

In order to complete the calculation of the Euclidean action, one has to compute the 
Einstein-Maxwell terms in (^]^). Since the scalar curvature TZ can be evaluated from the 
Einstein field equation 

^^^^u - l9^.un = 8^GT^, ^ 7^ = J^'^'J'^.u, (4.20) 

where 

V = 4;^(V-^- - 45m.-^°^-^«/3)> (4-21) 

is the energy-momentum tensor of the Maxwell field, it is enough to evaluate the electro- 
magnetic part of the action. 

For this purpose, recall that for a solution of the Maxwell equations J-'^^-u = 0, one 
can rewrite the Maxwell field Lagrangian density as J-^i,J-^^ = {2T^^ A^)-^u- Combining 
with ( |4.9| ) gives 

lib L = ^ / = 

We are now in a position to compute the Gibbs free energy. Substituting (|4T3| ), (|4lg| ), 
(|420D, and (|^) into yields 

-PG{mo,L) = -5^"" = 0. (4.23) 

Since G{mQ, L) = m — TS — $Q this implies 

m-^Q = TS, (4.24) 

where m is the total energy (mass) of the system. However, by the generalized Smarr's 
relation (d — 3)(m — ^Q) = {d — 2)TS + tL. Thus —TS = fL. Since all the quantities 
involved in the last equality are positive definite and L does not vanish, we conclude in the 
extremal case 

TS = f = 0, (4.25) 
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so the tension of the system is zero. Substituting this result back into ( [4.24| ) and using the 
extremahty condition (|3.9|), we see the mass of the system is unrenormahzed 



m 



= mo. 



(4.26) 



There remains one last thermodynamic quantity of interest. The redshift factor ( 3.14| ), 
is easily computed from our exact extremal solution and we have 



5oo(x) 



H-\x) R=^gooiO) = (1 + A)- 



(4.27) 



where in the last equality we excluded the n = term in the definition of H, since it 
corresponds to the short wavelength scale of order r^. This term is a pure infinity and 
is also present in the EFT approach, but in that case it is set to zero via dimensional 
regularization. For our purposes here, it can be ignored and thus does not contribute to 
the value of goo. Expanding the last equation in A, one recovers ( 3.18 ) for r/ = 0. 

The results in ( 4.25| ), ( 4.26 ), and ( 4.27 ) summarize the major outcomes of this section. 
In particular, these results verify the non-renormalization statements obtained perturba- 
tively in Section |3| via the EFT approach. 



5. Finite size effects 

In this section we add non-minimal couplings to account for finite size effects and study the 
resulting thermodynamics. We start with the uncharged case and then proceed to charged 
caged BHs. We determine the order the non-minimal operators enter the perturbative 
expansion for charged BHs. The existence of the exact solution in the extremal case allows 
us to undertake a matching calculation for the Wilson coefficients of the finite size operators 
in this limit. 



5.1 Uncharged caged BHs 

To incorporate the finite size operators into the action, we need to add all possible non- 
minimal operators which respect diffeomorphism and reparametrization invariance. These 
operators must use proper time throughout and can involve combinations of the geometric 
invariants, i.e. the Ricci scalar R, the Ricci tensor Rfj_u, and the Riemann tensor Rfj.upa- 
Higher derivative combinations of these terms are also possible. 

As usual in an EFT, the derivative expansion is truncated at a given order determined 
by the accuracy required in the calculation. In the uncharged case it was shown the 
operators involving R and Rfj^i^ can be removed by field redefinitions. These operators are 
redundant since the vacuum equations of motion are R^i, = 0. This means all physical 
finite size operators can be built from the Weyl tensor C^^p^ or alternatively from the 
Riemann tensor. Therefore, the leading order finite size operators are composed of the 
electric and magnetic components of the Weyl tensor squared [18]. Hence, the worldline 
action including these operators becomes 

SBu[x,g^,u] = -mo J dr + ji J drE^.E'"' + 6i j dTBp,,„p^_^B^''-^''^-\ (5.1) 
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Figure 6: Feynman diagrams with the finite size operator 71. (a) Contribution to the ADM mass, 
(b) Redshift modification, (c) and (d) Tension contributions from 71. 



where 



fix" dx^ 1 g dx°' dx^ 



Since our setup is static, the magnetic components of the Weyl tensor do not contribute 
to any thermodynamic observable and so we ignore the 61 operator. If we consider a non- 
static problem such as a scattering process, this operator may contribute since its Wilson 
coefficient 5i may be non-zero. 

The 71 operator will give rise to new vertices coupling to the worldline. The relevant 
worldline couplings will be (fxf) and (jxpOij. These are derived by expanding the operator 
and the result is 

^BH 13 -ii{didj(l)f - 2-iiaij{didk<p){djdk<t)), (5.3) 

where we have neglected terms which vanish by the leading order equations of motion. 
Terms with derivatives acting on Uzz in ( |5.3D have also been omitted since these do not 
contribute to our calculations. 

To derive the thermodynamics including finite size effects, we need to use a differ- 
ent approach compared to Section 0. There we assumed the local temperature Tq in 



( 3.11 ) is given by the temperature of a spherical uncompactified BH ( 3.12 ) and derived the 
asymptotic temperature using the redshift. This methodology is justified when the BH is 
spherical. However once finite size effects are allowed on the worldline, this is no longer 
the case. Instead, we must first calculate the mass and the tension of the system. With 
these relations we can derive the remaining thermodynamic quantities, independent of the 
BH shape. The redshift is also an observable so we compute it for completeness. 

A direct computation of the tension from Feynman diagrams requires we use the pre- 
scription of Q and compute J d'^^^xT^^ = —fL. We will however use the KK metric 
parametrization where the tension is obtained from diagrams with one external Gzz- In 
this case the sum of all diagrams yields — ^tLct^^. 

The Feynman diagrams needed to compute the contribution of the 71 operator to the 
mass, the redshift, and the tension are displayed in Fig. ^. Their evaluation requires some 
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new sum integrals which are hsted in Appendix We obtain 

Fig. |(a) = jimoid - l){d - 2){d - 3)V', (5.4) 

Fig. 1(6) = -271 (d -l){d- 2){d - 3)V', (5.5) 

Fig. 1(c) = -7imo(d -2){d- 3fid^ -id + 5)fi^a,„ (5.6) 

Fig. 1(d) = -27imo(d - 2f{d - ?>fi?a,,, (5.7) 

where we have defined 71 = 7i?TiorQ so 71 is dimensionless and the expansion parameter is 
defined by 



^1 d-i 



^=(T) C(d-l)~At^. (5.8) 



L 

The C'(^^) corrections to the mass, redshift and tension are then 

m = mo [1 - 7i((i - l){d -2){d- 3) V^l , (5.9) 
R = I - 2^i{d - l){d - 2){d - 3)V^ (5-10) 

— = 2^l{d-lf{d-2){d-2>f^?, (5.11) 
mo 

where we did not include the 0{\) and 0{}?) contributions from Section ^. We see that 
all corrections to thermodynamic quantities from the leading finite size effects scale as 
^2 ^ A , where ~ for d = 5, /x^ ~ for ti = 7 and /i^ ~ A^' with 2 < p < 3 for 
d>7. 

From the Smarr relation and the first law we find the entropy 

5 = 5o [1 + • 7i/i2] ^ (5_;L2) 

with Sq given by (|3.13| ) in the limit rj = \ and the temperature 

T = T{L ^ 00) [1 - 7i(3d - b){d -l){d- 2){d - 3)/x2] , (5.13) 

where T{L — > 00) is the uncompactified temperature of ( |3.12| ) with ry = 1. We note that 
when finite size effects are included, the entropy remains non-renormalized at leading order. 
Moreover, we observe that the redshift R does not coincide with the renormalization of the 
temperature - the local temperature To is not simply given by T(L co). In fact, we can 
now extract the local temperature including the leading order finite size operator to be 

To = T/R = T{L ^ 00) [1 - 71 (d + l){d - l){d - 2){d - 3,)^^] . (5.14) 

The Wilson coefficient 71 is still undetermined and its exact value must be extracted from 
a matching calculation at the scale tq. If possible, any matching calculation is best under- 
taken in d dimensions since there is no ambiguity when extracting the Wilson coefficients. 
In d dimensions the finite size contributions to the thermodynamics have a different di- 
mensional scaling compared to the point particle contributions. However, this is beyond 
the scope of the current work. Since 71 is dimensionless, one would expect 71 ~ 0(1); 
however, it can depend on the dimension d. 
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5.2 Including charge 

For a charged caged BH, we do not have a Ricci flat background, so we need to consider 
the finite size operators with two derivatives. In general, the effective BH action (| 
including all terms up to two derivatives is 



SB}i[Xfj,,g^,v,afj] = -mo j dr + Q j dx^a 



+ ai I drR + 02 / drR 



dx^' dx'' 



dr dr 



+ as y drf^^r + Oi,j drf^^f^ — — + ..., (5.15) 

where we have restored the time integrals for convenience and in the last two lines we have 
introduced the lowest order operators describing finite size effects. According to the power 
counting rules ||8|, [l^, the Wilson coefficients of these operators scale as mor^ if we power 
count Q ~ tuq. 

Due to the presence of electromagnetism, operators involving the Ricci tensor with 
coefficients ai and Q2 do not vanish by the leading order equations of motion. However 
we can still use the equations of motion, where i? ~ // for the Ricci tensor and scalar, see 



( 4.20 ). Thus we can eliminate the ai and 02 operators in favor of the 03 and 0:4 operators. 
Alternatively, it is also possible to redefine the coordinate system to eliminate 03 and 04, 
and the procedure is similar to the uncharged case discussed in j^. 

We will eliminate the ai and 02 operators and use a different operator basis for the 



electromagnetic finite size operators than in ( 5.15 ) where the and operators are 
expressed in terms of the electric and the magnetic b^^ ^^^^^_^ components of the Maxwell 
field strength tensor 



dx" 

d7 ' 

1 dx"^ 

^Mi---Md-3 ~ (rf^r3yT~^^"'^i---^'*-3^'>''^^''^ ■ (5.16) 

For a static system the operator formed from the magnetic field will not contribute to the 
thermodynamics, so we will neglect it. The worldline action including the leading order 
static finite size operator with two derivatives reads 

SBYL[x^i,g^_iu,a^] = -mo j dr + Q j dx^'a^ + a j dre^e^ + . . . . (5.17) 

We now compute the effect of the a operator. In the BH's rest frame the a operator 
reduces to —e'^'^^^'^~^^j^^{di(p){dj(p). The Feynman rule for the vertex coupling two (/^'s to 
the worldline is 2ak • p, with both momenta incoming. The corresponding mass renor- 
malization diagram is shown in Fig. 0(a), which scales as ~ moA2('='-2)/('^-3). Calculating 
this diagram shows that it vanishes, and this has a clear interpretation. The finite size 
operator in ( ^.17 ) encodes deformation due to the electric dipole polarization of the BH. If 



we think of the problem as a line of charged BHs at regular intervals in the covering space. 
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Figure 7: (a) Contribution to the ADM mass m from the finite size operator insertion. The open 
square denotes the vertex resulting from a J dre^ic'' insertion, (b) and (c) Tension contributions 
from a operator, (d) and (e) ADM mass renormalization diagrams resulting from the quadrupole 
finite size operators. The solid square vertex denotes the appropriate quadrupole finite size opera- 
tors, (f) Contribution to the rcdshift of operator 72. (g)-(k) Tension contributions from the 72 and 
73 finite size operators. 



as in Fig. H, then the polarizing effect from a charge to the right will be canceled by the 
equidistant charge on the left. Therefore the overall polarizing effect is zero and there is 
no mass renormalization from this worldline operator. Similarly, the diagrams involving 
this operator and an external azz which contribute to the tension, see Fig. ^(b) and (c), 
are found to vanish. 

As a result, the first non- vanishing finite size contributions to the BH thermodynamics 
come from operators with four derivatives, which encode deformations due to induced 
quadrupole moments. As explained above, in a static system we can form all relevant 
operators using the electric field e^, the electric components of the Weyl tensor E^^, and 
higher derivative combinations. The static finite size terms in the action read at order four 
derivatives 

.5bh[x, a^] =71 J drE^.E^''' -j2 J (ir(V^e,)^^^ + 73 J (ir(V^e,)(V^e'^), (5.18) 
where the first operator coincides with the one in the uncharged case of (f 



The finite size terms in the action (5.1^) lead to new worldline vertices with two scalar 



fields, either gravitational or electromagnetic. There are also vertices on the worldline with 
two scalars and one azz which contributes to the tension calculation. The leading order 



correction to the ADM mass due to an insertion of these operators introduced in ( 5.1 



arises from the diagrams Fig. ^(a) and Fig. ^(d) and (e). Calculating the new diagrams 
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gives the renormalized mass 
m = niQ 1 



1 Q d-2 l/Q\V(i-2^^ 
71+72- ; — ^+73 



2 rriQ d — 3 



4 V^o / \d — 3 



{d-l){d-2){d-3fij.'^j , 
(5.19) 



where we defined the dimensionless couplings 7^ from 7j = ^imorQ. 

For the redshift factor R we compute Fig. ^(6) and Fig. |^(/) with the result 

= 1 - ( 271 + j^l — ^) id -l){d- 2){d - 3ff?, 



(5.20) 



' 2 rriQ d — 3 

and the tension is computed via the diagrams Fig. ^(c) and (d) and Fig. 0(5')-(A;) to be 
fL 



mo 



1 Q d-2 l/Q\Vd-2^^ 

71 +727T J ^+73 

2 mo d — 3 



4 \mo / \d — 3 



{d-lY{d-2){d-3Yfi'. (5.21) 



In the charged case we have two dimensionful short distance quantities, mg and Q, so in 
general the dimensionless couplings 7^ are functions of the order parameter ij. Thus the 
thermodynamics derivation of the remaining quantities S, T and $ becomes cumbersome 
and not insightful since we have to keep derivatives of the couplings with respect to Q and 
mo. 

In the extremal limit = 0, we can use the exact solution of Section |^ to undertake a 
matching calculation for the undetermined Wilson coefficients. This is done by comparing 
the exact results m = mg, f = 0, and R = (l + X)-^ to the expressions in the perturbative 
expansions ( 5.19| ), ( |5.2C| ), and ( p.21| ). The non-renormalization of the mass and tension are 
seen to yield one relation between the Wilson coefficients 7^, since the linear combination 
of coefficients in (|5.19|) and ( ^.21 ) is degenerate. 

The remaining relation is derived from R. We note the 0{ijl^) contribution must vanish, 
since it is non-analytic in A in arbitrary dimensions. We find from ( 5.2C| ) 



72 

71 



r;=0 



,d-3 
'd-2 



sign(Q). 



(5.22) 



Using this relation and m = mo or f = we get 



73 

71 



?7=0 



d-3 
'd-2' 



(5.23) 



This constitutes a partial matching of the Wilson coefficients of the finite size operators in 
the extremal limit. 

The matching calculation presented above has several shortcomings. The calculation 
of the thermodynamic properties of caged BHs cannot yield the coefficient of any operator 
involving magnetic type components, since these do not contribute for static systems. 
Further, even in the case of the electric type operators, the leading finite size operator with 
two derivatives e^e^ does not contribute to the thermodynamics of caged BHs - so we 
cannot fix its Wilson coefficient a. Lastly, while the three electric type operators with four 
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derivatives do contribute to the thermodynamics, we are not able to fix all three Wilson 
coefficients 7^ uniquely, due to the degeneracy between the mass and tension relationships. 

We have not exploited the non-renormalization of the electrostatic potential We 
might be able to fix the Wilson coefficients 7^ uniquely if we could calculate the electrostatic 
potential $ in the EFT. If we succeeded in such a computation, the 0{fj,'^) contributions 
from the finite size operators would most likely entail a linear combination of 72 and 73. If 
this new equation was linearly independent of the other equations derived here, this would 
imply 71 = 72 = 73 = in the extremal limit. 

6. Conclusion 

Using an EFT approach we have analyzed the thermodynamic properties of small com- 
pactified BHs carrying charge. We obtain the relevant thermodynamic quantities to O(A^). 
Standard power counting arguments show that up to O(A^), all thermodynamic contribu- 
tions arise from a point particle description of the BH. This implies that to this order the 
horizon is spherical and not deformed due to the presence of the compact dimension, and 
the entropy is given by the area of a spherical horizon as for an uncompactified charged 
BH. 

In the extremal limit, we constructed an exact solution using the standard methods of 
GR. From the exact solution, we find that the mass m, tension r, and electrostatic potential 
$ are non-renormalized in the extremal limit, and the redshift is given by a geometric series 
in A. We use the extremal thermodynamics and the uncharged perturbative results of [^, |l2| 
as checks of our perturbative results for the thermodynamic properties and find agreement. 

The leading finite size corrections to thermodynamic properties are computed in both 
the charged and uncharged cases in terms of Wilson coefficients of higher order operators. 
These coefficients need to be determined by a matching calculation. In the uncharged case 
the entropy does not acquire a correction due to the leading order finite size operators. 
If the Wilson coefficient is non-zero, this is interesting, since this means its area remains 
the same even though the horizon is now deformed and no longer spherical. Moreover, the 
local temperature and the temperature of an uncompactified BH do not coincide when we 
include finite size effects, which is a clear sign of a non-sphericity. 

For charged caged BHs, we find the leading finite size operators with two derivatives 
do not contribute to the thermodynamics. Also, due to the presence of two short distance 
scales Q and niQ, the dimensionless Wilson coefficients can plausibly depend on their ra- 
tio. We computed the leading finite size corrections to m, R, and f, which come from 
four derivative operators. In the extremal limit, we undertake a matching calculation by 
comparing with the exact solution and we constrain two out of the three Wilson coefficients. 

Our results contribute to the study of the phase diagram for black objects with one 
compactified dimension. In particular, for dimensions d > 7, the leading finite size opera- 
tors yield the dominant correction to the existing perturbative results at 0{X'^). It would 
be interesting to perform a complete matching of the Wilson coefficients of the finite size 
operators. For the uncharged case, it may in fact be possible to obtain an estimate for 71 
from numerical data for d > 7. 



- 21 - 



In the uncharged case, Q performed an interesting comparison between their pertur- 
bative results and the numerical results of Q. This comparison shows the perturbative 
result at 0{X) matches the numerical solution to 10-20% all the way up to the point where 
the BH/non- uniform BS phase transition occurs. Since we have shown the finite size effects 
first enter at the same order as in the uncharged case, we can expect our leading order 
results to behave similarly. We can only make general inferences at this point, since the 
details of a possible transition are not captured by the EFT. 

Based on our result for the leading order tension f oc r/^A, for a given bare mass mo and 
compactification L, the addition of charge to the caged BH will reduce the tension. Due 
to the excellent agreement between the EFT and numerical treatments of the uncharged 
caged BH, we speculate that this property will also hold non-perturbatively. Therefore one 
can reasonably expect the topology changing transition to a BS phase to be delayed when 
charge is added to the BH. 

Intuitively this increased stability of the BH phase can be understood in the covering 
space. Here the transition to a BS will occur when we increase the BH mass and the 
horizons of neighboring image BHs overlap. If there was no horizon deformation this 
would occur at values of the bare mass when the horizon radius is r = L/2. However, the 
attraction between the neighboring image BHs can cause this overlap to occur at lower 
bare mass due to horizon deformation. Since like charges repel, the tension from the 
gravitational attraction will be reduced by an electrostatic repulsion of the image BHs - 
making a charged caged BH relatively more stable than an uncharged caged BH with the 
same bare mass. 

So far there have been few studies of the BH/BS phase transition with the inclusion 
of charge. We have covered the entire mass-charge phase space for small A and our study 
provides a first exploration of the BH side of the phase space. Due to the presence of the 
additional parameter Q, a richer structure of the phase diagram can be expected. Further 
study, either numerically or analytically, towards completion of this phase diagram, is left 
for future work. Investigation of the possible phases in the extremal or near extremal limit 
may exhibit new phenomenology. 
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A. Table of Feynman integral sums 

All integrals are derived using the techniques described in |3^], with the sums being Rie- 
mann zeta functions. The first two results listed below, Io{L) and Ii{L), were derived in 
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Q and I2{L) was computed in [12|. We list all integral sums used in this work below 
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B. Thermodynamics via the Helmholtz free energy 



Here we demonstrate the calculation of the thermodynamics of Section using the 
Helmholtz free energy F, without using any assumptions about the entropy S of the charged 
caged BH. We start by defining the Helmholtz free energy 



m — TS. 



(B.l) 



The first law (|3l2ll ) gives the differential relation dF = —SdT + '^dQ + fdL. We now derive 
the following differential relations between the Helmholtz free energy and the entropy S, 
the electromagnetic potential <I>, and the tension f, 



^ " 'or 



$ = + 



r = + 



+ 



f8F\ 


/ dT 






d{F,L,T) 


' f dT 


d{mo,L,Q) 




d{F,Q,T) 


Y dT 


d{mo,L,Q) 





0/ L,C 



(B.2) 
(B.3) 
(B.4) 



After the second equality in the above relations, we have rewritten the thermodynamic 
relations to be functions of mo, L, and Q, since these are the bare parameters of the EFT. 

At this stage we use the Smarr relation ( ^.24 ) to rewrite F eliminating the temperature- 
entropy term. This gives 



{d-2)F = m+ {d- 3)^Q + fL. 



(B.5) 
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This form is now amendable to solution, since we can form a differential relation for F 
through the above equations for <I> and r, ( |B.3| ) and ( |B.^ ) respectively. Doing this allows 
us to rewrite ( [B.5| ) as follows 



{d-2)F = m+ [{?,- d)Q 



d{F,L,T) 



+ L 



d{F,Q,T) 



d{mQ,L,Q) d{mo,L,Q) 



dT 



dmr 



0/ L,i 



(B.6) 



We have now obtained what is obviously a first order multi-variable partial differential 
equation for the Helmholtz free energy F. A solution can be constructed for F through 
iteration, first order by order in A and then order by order in 77 at a given order in A. After 
some calculation, we find F to 0{X'^), and this is given by. 



F = niQ 



d-3 
d-2 



V + 



d-3 
d-2 



1 



7] 



d-3 



1 



d 



— 7? H 7?" 

2 ' 2 ' 



ld-3 



2d-2 



A^ 



(B.7) 

The entropy S, the electrostatic potential ^, and the tension f are computed using the 
thermodynamic relations in ( [B.2D -( |r!4 ), and their results agree with the ones stated in 
Section 3.2. Using a Legendre transformation, the Gibbs free energy is G = F — and 
we find agreement with ( p. 27 ). 
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